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The effects of elastic foundation on the fundamental 
natural frequencies of isotropic plates and orthotropic 
laminates are investigated, A finite element model is deve- 
loped to take into account the shear rotation and rotary 
inertia effects. A generalised foundation model based on 
the concept of isotropic elastic half space is used. This 
model can be reduced to one parameter or two parameter 
models. Simply supported square plate/laminate is considered 
■ for the analysis. Effects of varying foundation stiffness, 
ply angles, aspect ratio and side to thickness ratios on the 
fundamental frequencies are investigated. The results are 
plotted and tabulated. 



CHAPTER 1 
INTRODUCTIOM 

Rectangular plate is one of the most commonly used 
structural members in present day's industrialised world. 

In most design problems ensuring the safety of the plate by 
mere static analysis is often inadequate. Therefore, it 
requires a dynamic analysis that takes into account the tim.e 
varying forces, periodic t)r random, as experienced in cases 
of rectangular plates used on ships and submarines, subjected 
to a tangential liquid force or those experienced by the off- 
shore drilling platforms due to fluid plate interaction. 

Since the advent of fiber composites, composite 
plates and laminates are fast replacing those made of con- 
ventional metals and alloys in specialised weight sensitive 
applications like aerospace and marine. In general, composite 
materials are most ideally suited and are cost effective for 
structural applications where high strength to weight and 
high stiffness to weight ratios are required. In the recent 
years, development of solid propellent rocket motors, use 
of soft filaments in aerospace applications intensified the 
need for solutions of various problems of beams, plates and 
shells continuously supported by elastic or visco elastic 
media. 

The fact that the free vibration analysis for funda- 
mental natural frequencies is the first logical step for 
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dynamic analysis of a structure, need not be overemphasized. 

In the present study fundamental frequencies of a composite 
laminate resting on an elastic foundation are investigated. 

1.1 Literature Review: 

Since Kirchoff's thin plate theory, which neglects 
the effects of rotary inertia and transverse shear, a number 
of successful attempts have been made to improve the theory by 
introducing the above effects. One of the most widely accep- 
ted among these is Reissner's work [l j which included the 
effect of transverse shear deformation in the plate bending 
theory. Later, R.D. Mindlin [2] included the effect of rotary 
inertia in addition to the transverse shear effect in flexural 
motion of isotropic plates. Analytical solutions for static 
analysis of thick plates with regular geometries and boundary 
conditions have since been obtained by various authors, in 
particular D.H. Donnel et al [3], Salerno and Goldberg [4] and 
Carley and Langhaar [s]. Dynamic analyses of thick rectan- 
gular plates were carried out by Mindlin, Shacknow and 
Deresiewicz [6] which included the effects of rotary inertia 
and transverse shear in flexural vibrations. Srinivas and 
Rao [ 7 ] presented an analysis of bending, buckling and vibra- 
tion of simply supported thick orthotropic rectangular plates 
and laminates by classical theory and compares the exact 
results with those obtained by thin plate theory in cases of 
displacements, stresses and natural frequencies. Smith [s] 
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analysed moderately thick rectangular plates for bending, 
using rectangular finite elements with nine generalised 
displacements at each node. Grieman and Lynn [9] developed 
a rectangular finite element model with transverse displace- 
ment and two orthogonal rotations as the nodal degrees of 
freedom. Pryor, Baker and Frederick [lOj formulated finite 
element bending analysis of Reissner plates using rectangular 
elements with four nodes, each node having five nodal degrees 
of freedom viz., transverse displacement, two transverse 
shear strains and two total rotations. 

In the field of finite element analysis of plate 
dynamics, a significant contribution is made by Rock and 
Hinton [11, 12], who introduced an 8-noded rectangular 
isoparametric element to take into account the transverse 
shear deformation. They also suggested different schemes 
for lumping masses to obtain the mass matrices that take into 
account the effect of rotary inertia. Reddy [13] formulated 
the free vibration analysis of laminated plates by finite 
element method using the theory postulated by Yang, Norris 
and Stavsky [14]. 

A critical study of various foundation models was 
presented by Kerr [is]. A finite element method to incorpo- 
rate the foundation stiffness in the structure stiffness was 
first suggested by Cheung and Zienkiewicz [16], in which 
static analysis of plates and tanks resting on elastic foun- 
dations was carried out. This work was later followed up by 
Cheung and Nag [17] who modelled the elastic continuum of 
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the foundation as isotropic elastic half plane and elastic 
half space to analyse the linear and non linear behaviour of 
beams and plates on elastic foundations. Svec [18] investi- 
gated thick plates on Winkler foundation. Kameswara Rao 
et al investigated dynamic response of beams on generalised 
elastic foundations [19]. More recently, Thangambabu et al 
[20] carried out frequency analysis of thick orthotropic 
plates on Winkler foundation by finite element method using 
high precision triangular elements. 

1.2 Objective and Scope of Present Works 

To date not much work has been published on the dyna- 
mic analysis of plates resting on elastic foundations. In 
the present analysis, an effort is made to investigate the 
effects of elastic foundation on the fundamental natural 
frequencies of laminates by finite element method. The study 
is limited to small amplitude vibrations, preserving the 
linearty. Effects of varying ply angle, aspect ratio, side 
to thickness ratio and foundation modulus on the fundamental 
frequencies are tabulated. 
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CHAPTER 2 

ANALYTICAL FOFyAULATIW OF THE PROBLEM 

2.1 Introduction: 

In this chapter, the concepts of anisotropy and 
orthotropy and their application to derive the laminate 
stress strain relationships are briefly introduced prior to 
the analytical formulation of laminate stress strain rela- 
tionships. Effect of shear deformation is included in the 
laminate stiffness and the resultant equations of motion 
are presented. Formulation of the elastic foundation is 
presented at the end. 

2.2 Stress Strain Relationship for Anisotropic Materials: 

The generalised Hooke's law relating the stresses 
and strains for an anisotropic material with no planes of 
symmetry for material properties is given by 
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By introducing corresponding planes of symmetry of 
material properties in (2.1) these relationships for 
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different cases of anisotropy can be obtained. The stress 
strain relationships for an orthotropic material in prin- 
cipal material directions can be obtained as. 
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which has nine independent elastic constants, free from 
interactions. 

These elastic properties are defined in terms of 
engineering constants E|^, E^, G^t principal material 

directions [21]. More general form of these relationships 
for anisotropic bodies can be found in reference [22]. The 
relationships defined in (2.2) can be transformed to any 
arbitrary orientation of fiber in case of continuous fiber 
lamina. 


{o) = [q] Ce) (2.3) 

where [Q] is the transformed stiffness matrix. 

2.3 Application to a Laminates 

2.3.1 Laminate strain displacement relationships: 

A laminate is two or more laminae in arbitrary fiber 
orientation bonded together to act as an integral structural 
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element. The laminate is presumed to consist of perfectly 
bonded laminae with non shear deformable bonds, so that the 
displacements are continuous across the contact area. We 
assume that the normal to the middle surface remains straight, 
unaltered in length (i.e. = 0) and normal under deformation 

The displacements u and v at any point z through 
laminate thickness is (Figure 2.1) 


u 


= u_ - z 


= v_ - z 


9w 
o 

3x 

3w 

g 

3y 


(2.4) 


For small strains, the strain displacement relations are 
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Similarly the mid plane strains in terms of displa- 
cements can be defined as 
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Using (2.4) through (2.6) the total strains can be expressed 
in terms of midplane strains and curvatures as 
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2.3.2 Resultant laminate forces and moments: 

The resultant forces and moments per unit length 
acting on the laminate are obtained by integration of the 
stresses in each layer/lamina through the laminate thickness 
For an N layered laminate v/ith thickness h, the stress and 
moment resultants are obtained as 
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and 



( 2 . 10 ) 

Using the transformed stiffness matrix of equation 
(2.3) in the above, the force and moment relationships are 
obtained as 
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2.4 Formulation of Laminate Stiffness with Shear Deformation: 


2.4.1 Strain displacement relations: 

A laminate with constant thickness composed of even 
number of thin anisotropic layers oriented at 8 and -0 alter- 
natively, is considered with the origin of the coordinate 
system located at the middle plane and z axis being normal to 
the middle plane. 

By Yang, Norris and Stavskey’ s theory [l4] the 
assumed displacement field in terms of midplane deformations 
and shear rotations 

u = u^(x, y, t) + z . WyCx, y, t) 

o ^ 

V = v^(x, y, t) + z . ^„(x, y, t) (2.13) 

o y 

W = w(x, y, t) 

Applying the strain displacement relationships from 
equation (2.5) to include the shear rotations, we get. 
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2.4.2 Resultant laminate stiffness: 

The constitutive relationships for an orthotropic 
layer in principal material directions L and T are given by 
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Now, using the relationships for resultant forces and moments 
for a laminate from equations (2.4) through (2.10), we can 
express the constitutive equation (2.14) in terms of forces 



We can assume that the shear rotation terms are independent 
of each other and make = 0, i.e. the coupling between 
and Qy is neglected. 


2.4.3 Equations of motion: 

The equations of motion are obtained by v/riting the 
expressions for the kinetic energy T and potential energy V 
and taking the first variation of the Lagrangian L = T - V, 


as 
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2.5 Formulation of Elastic Foundation: 


2.5.1 Foundation flexibility and stiffness matrices: 

Boussinesq [23] formulated the deflection formulae 
for isotropic elastic half space for point loads. Treating 
the forces acting at nodes as equivalent uniform pressure 
acting on a rectangle around the node, t..e coefficients of 
the foundation flexibility matrix are obtained [l7] as. 


for vertical pressure (Figure 2.4): 


V„{1 - yf ) 1 . 1 

„ _ JG — L_ (b sinh” g- + sinh^'B - C sinh 

mn btiEj, d 


^ ^ - sinh"^C) 




FIG. 2. 4 VERTICAL 
OF ELAS' 




18 


u. 


mn 


V^(1 - 2a'^)(l + 


2ti E, 


[B tan"'' i - In — — ^ 

4T? 


and 


G tan ^ i + In — ] , m n 


vT- 2 
1 + C 


w 


mn 


2V 

[B sinh"'' g- + sinh“^ B](1 - 


u =0 
mn 


m = n 


for horizontal pressure; 


H (1 - 


u 


n 




— [B sinh”'' K + sinh“''B - C sinh"'' 4 - sinh ^ C J 


a B 


H 3/ (1 + 3^ ) 

+ -2 — — - ■ ■ - - [B sinh“' B - sinh” C] , 


m n 


2H (1 - yt) 


u 


n 


MN 


u E^ a 


— [B sinh"^ ^ + sinh"^B] 


2H y. (1 + 3^.) . 

^ — Q— I— L- sinh-'B 

-K E^ a 


w, 




u 


mn 


m = n 
(2.19) 


where, 

B = [2(iii - n) + 1 ] f 
C = [2 {ib - n) - l] f 

a, b = sides of rectangle around the node. 

The coefficients thus obtained from the equations 
(2.19) form the flexibility matrix of the foundation. This 
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matrix is inverted to obtain the stiffness matrix of the 
foundation. 
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2.5.2 Generalised foundation stiffness: 

The stiffness matrix obtained in equation (2.20) is 
with respect to the foundation plate interface and that of 
the plate in general, is with respect to the neutral surface. 
Therefore the terms involving the horizontal pressures need 
to be modified to account for the plate thickness (Figure 2.5 ). 

At any node i, the relationships between plate and 
foundation displacements are given by 
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with primed symbols representing foundation system. 

Using (2.21 ) in (2.20) and expanding, we obtain the 
generalised stiffness matrix which takes into account the 
moment rotation terms as follows. 
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This equation 

for 

generalised 

foundation 

model 

can 


be modified to obtain either a one parameter model (Winkler 
foundation) or a two parameter model with interactions bet- 
ween vertical displacements by supressing the unnecessary 
terms. 


2.6 Overall Stiffness Matrix; 

The overall stiffness matrix for the laminate resting 
on elastic foundation is obtained by adding the foundation 
stiffness matrix to that of the laminate. 
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CHAPTER 3 

FINITE ELEV^ENT FOR^.'.ULATICW 

3.1 Introduction: 

In general plate bending analysis in the case of 
isotropic plates is done by considering the transverse dis- 
placement and the tvjo shear rotations as field variables. 

In the case of fiber composites, the transverse shear effects 
are more pronounced in comparison with the isotropic case as 
the transverse shear modulus is very small compared to the 
inplane Young's modulus. The inplane displacem.ents have a 
significant effect on the natural frequencies of composite 
laminates. In this chapter a finite element model for the 
dynamic analysis of composite laminates is developed. The 
effect of elastic foundation is incorporated in the plate 
stiffness matrix. 

3.2 The Element: 

An isoparametric quadrilateral plate bending element 

with u , V , w, ? , f as the field variables over the element 
o o X y 

is chosen for the analysis (Figure 3.1). 

3.2.1 Shape functions: 

Shape functions associated with each node define, 

(i) geometry of the element in terms of nodal coordinates 
(ii) displacement variation over the element in terms of 
nodal displacements. 
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At any point within the elei-ient 
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The shape functions associated with each node for the 
element in terms of non dimensional! sed coordinates r/ are 
as follows. 
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Same shape functions are employed for all the five 
nodal displacements. 

3.3 Elemental Matrices: 


3.3.1 Strain displacement relationships: 

The generalised strain vector in terms of displacements' 
from the constitutive equation (2.17) is 
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where 


[B] = [L] [N] 


(3.6) 


L is the linear operator matrix obtained from the nodal strain 
vector (3.4) as 
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3.3.2 Stiffness matrix: 

Stiffness matrix is given by [24] 
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^ ~ [D] [B] drj det J 


(3.8) 


with det J accounting for the transformation of the cooiciinates 
from local to global and J is defined by 


[J] = 


ox 

a1 

9x 

3h 


ti 

3n 

iY 

oh 


(3.9) 
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(3.10) 


[D] is the material property matrix occurring in the stress 
strain relationships. 

The numerical integration in equation (3.8) is carried 
out by 2 X 2 Gauss quadrature [1l]. 


3.3.3 Mass matrix: 

Consistent mass matrix is derived using the same 
interpolation functions as given in equation (3.3). 

Nodal accelerations at any node r are given by 





(3.11) 


u V and w produce inplane and lateral inertia forces re spec- 
0^0 

tively and i , ^ produce rotary inertia couples. The 

y 

inertia force vector can be written as 
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{a)^ = {P (I P V , p w, ^ f ^ 

O’ o* ^ X_| 


(3.12) 


12 X* 12 y 

At any point in the element, inertia force can be exoressed as 


whe re 


e it 


V6^ } 


fa 1 = [p] [N ] {6® ] 


^2 * • • • '^ 8 ^ 


and 


0 


0 


sym 


Ph^ 

12 


0 0 

[ P ] = Poo 

Oh^ 

0 

Ph^ 

12 J 

Nodal forces due to inertia can be expressed as 


F? = / [N] {a } dA 


Using (3.13) in (3.15) 

= / [N]^ [p] [n] dA ! 

= [hf] fS®) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 

(3.17) 


where 

[M® ] is the elemental consistent mass matrix given by 
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[M®] = / / [n]^ [P] [n] dS dri det J (3.18) 

-1 -1 

The consistent mass matrix has been evaluated using 
3x3 Gauss quadrature rule [11 ]. 

3.3.4 Resultant stiffness matrix: 

The resultant elemental stiffness matrix is obtained 
by adding the foundation stiffness of equation (2.20) to the 
elemental stiffness obtained in equation (3.8). 

3.4 Assembly of Elemental Matrices: 

Assembly of elemental matrices into global matrices 
is done by imposing the conditions that the displacement 
vectors are common and sum of the generalised force vectors 
are zero at common nodes. 

3.5 Solution for Free Response: 

With no forcing and damping terms existing, the 
dynamic problem is defined as 

[m] tsi + [K] (A! = 0 (3.19) 

By substituting tA)= (Z le^ in (3.19) we obtain the typical 
eigen value problem 

([k] - [m])U> = 0 (3.20) 

For the non zero solutions of the above eigen problem, the 
determinant of the equation (3.20) must be zero. 
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I [K] _ [?/j { = 0 (3.21 } 

This ecjuation yields the n values of v.'hen the size of 
matrices [K] and [m] is n x n. Usually, in structural problems 

[K] and [M] are positive definite thereby resulting in positive 
eigen values. 



4.1 Introduction 
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. f ■. ' *~5, yi, 
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A computer program based on the finite element formul- 
ation presented in Chapter 3 is developed to investigate the 
effects of elastic foundation on the fundamental natural frequ- 
encies of isotropic plates and orthotropic laminates. Two 
types of orthotropic plates ere used, biaxial s'/rumetry of the 
plate is made use of and quarter plate model is used to compute 
the fundamental frequencies and a half plate model is used to 
compute frequencies at higher mooes. The mesh is of the 
order 2x2 for the quarter plate and that for the half plate 
is 2 X 4. Simply supported boundary conditions are shorn in 

Figure 4.1. 


4.1.1 Material properties: 

Following material properties are used. 

Isotropic case: 

. , „6 

H = 30 X 10 psi 

j; = 0.3 

p = 0.283 Ib/in- 

Orthotropic case (graphite epoxy): 

6 

E = 27 X 1 0 psi 
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Ej = 0.675 X 10^ psi 

Glt = 0.405 X 10^ psi 

^'LT = 0.25 

P = 0.056 Ib/in^ 

Four layers of equal thickness made of same orthotropic 
material is assumed to constitute the laminate. 

Calculation of D matrix is based on the equation 
(2.17) for orthotropic case and the standard stress strain 
relationships are used for the isotropic case. 

4.1.2 Accuracy of the model: 

Numerical accuracy of the program developed is estab- 
lished by comparing the results obtained in earlier works 
6, 11, 17 . Comparisons in respect of follovdng are made. 

i) fundamental frequencies of a simply supported, 
square isotropic plate (Table 1) 

ii ) fundamental frequencies of square simply supported 
ortho tropic laminate (Table 2) 

iii) static deflection of central point of square 
isotropic plate on “Winkle r foundation (Table 3). 

In the third case, the generalised foundation model is 
reduced to that of a Winkler type and the results are compared. 
All tabulated frequencies are non dimensionalised as 

shovm . 

4.2 Isotropic Plate: 

One parameter, two parameter and generalised foundation 
models are used to investigate the effects of elastic foundation 


32 






I 


Table 1. 


Isotropic non dimensional 

,2 1/2 




H - ^ » h 


frequencies 

=10, simply supported 


Mode 

Mindlin [6] 

Rock & Hinton [l2] 

Present model 

1. 1 

0.0930 

0.0931 

0.09313 

2, 1 

0.2217 

0.2237 

0.2239 

3, 1 

0.4144 

0.4312 

0.4317 

3, 2 

0.5197 

0.5379 

0.5500 
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Table 2. 


Orthotropic non dimensional frequencies D 
simply supported plate 


_2, P 

wa ( , 

p h 


Q 

1 — 

1 

t 

-4 

300 

■T 

! 

! 

1 

450 

a/b 

a/h 

I 

I 1.0 

1 

t 

1 .2 

“1 — 

f 

i 1.0 

t 

! 

1 

1.2 

10 

17.616 

*(17.689) 

19.810 

(19.880) 

18.360 

(18.609) 

21 .495 
(21.567) 

20 

21 . 344 
(21 .281 ) 

24.427 

(24.362) 

22.385 

(22.584) 

26.877 

(26.837) 


* Ref. [13]. 


Table 3. Central point displacements of simply supported square 
isotrooic plate on Winkler foundation 



1 

h/a 1 




! 0.03 

0.05 

0.10 

E/E^ 

f 

1 

f 

1 ...... 



10 

0.56 X 10"^ 

0.41 X 10"^ 

0.28 X 10"^ 

*(0.79 X 10"^) 

(0.49 X 10"^) 

(0.25 X 10"^) 


0.37 X 10"^ 

0.218 X 10“^ 

0.1025 X 10"^ 

100 

♦(0.325 X 10“^) 

(0.22 X 10“^) 

(0.120 X 10”^) 


* Ref. [is]. 
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on the fundamental natural frequencies. The results obtained 

by using the two parameter and generalised foundation models 

are found to almost coincide. 

The fundamental natural frequencies are found to increase 

E 

with the increase in foundation stiffness. At values of 

tf 

smaller than 100 (very stiff foundation) the fundamental frequ- 
ency increases sharply as shown in Figure 4.1 

The mode shapes are investigated and it is observed 
that all the modes are predominantly bending modes. The 

results obtained by using the two parameter foundation model 
are lower than those obtained by the one parameter model. 


4.3 Symmetric Laminate: 

A symmetric orthotropic laminate of the type is 

used to investigate the effects of varying ply angle and the 
ratio EL/Ef fundamental natural frequencies of a simply 

supported case. For relatively soft foundations the fundamental 
frequencies are found to increase with increasing ply angle from 
Oo to 45 and are symmetric with respect to uhe 45'^ ply» 
frequencies increase with decreasing E^/E^ ratio. The bending 
modes dominate in softer foundations as well as in the stiffer 
foundations <100). However at E^/E^ ratio greater than 50. 

the fundamental frequency at 45- ply angle Is found to have a 
smaller natural frequency than those at other ply orientations. 
The smallest frequencies do not vary with further Increase in 











oo 

sj> 


the foundation stiffness. The results are plotted in Figures 
4.2 and 4.3. 

4.4 Antisymmetric Laminate: 

An antisymmetric angle ply laminate of configuration 
[+0/-0/ + 0/ - 0]^ is chosen. As in the previous cases, 
investigations were carried out for a square simply supported 
laminate. 

At lower foundation stiffness (higher values of 
the natural frequencies vary as in the case of the symmetric 
laminate. The predominant mode is the bending mode. With 
stiffer foundations the stretching mode is excited due to higher 
bending stiffnesses. However, a sharp variation in funda,;.enta_ 
frequencies is observed around E/E^ ratio of 10. The highest 
frequency belonging to the 30^ ply orientation and the smallest 
belonging to the 45 ^ ply orientation. These frequencies correspond 
to the stretching modes. The smallest frequency in the stret- 
ching mode at 45 c orientation can be attributed to the smallest 
values of E^ and (in plane moduli) at 45° ply orientation. 
However, this behaviour is restricted to this type of laminates 
only and separate investigations are necessary for difieient 
types of laminates as the laminate stiffness ixseli is affected 
by the change of any of the laminate parameters. The behaviour 

is plotted in Figure 4.4. 

The parametric effects are tabulated by varying the 
aspect raUo and side to thickness ratio at different ply 
angles and the results are tabulated in Table 4. The van-^io 






Table 4. Parametric effects: simply supported antisymmetric angle ply laminate 
(2 parameter foundation) 
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of fundamental frec^encies at different ply angles and paramo- 
trie ratios are plotted in Figure 4.5, 

4.5 Conclusions: 

Since the laminates can be made to suit the require- 
of desired anisotropy, separate analysis wherever is 
required is to be carried out to ascertain the dynamic beha- 
viour of the particular type of laminate. An extension of 
p-esent work can be undertaken for non linear vibration 
analysis including damping and for other boundary and external 
loading conditions. Experimental work to ascertain the beha- 
viour of laminates on elastic foundations may be undertaken. 

It is hoped that the present work can form a basis for under- 
taking the dynamic analysis of plates that can incorporate the 
interaction of the inertia of foundation system with that of 
the structure. 
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